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The interaction cross sections for ^Na + 12 C reaction are calculated 

using Glauber model. The continuum Hartree-Bogoliubov theory has been 



generalized to treat the odd particle system and take the continuum into 
account. The theory reproduces the experimental result quite well. The 
matter distributions from the proton drip line to the neutron drip line in 

>. 

^sO Na isotopes have been systematically studied and presented. The relation 

between the shell effects and the halos has been examined. The tail of the 

qq | matter distribution shows a strong dependence on the shell structure. The 

Os" 

neutron N = 28 closure shell fails to appear due to the coming down of the 
2p 3 / 2 aR d 2pi/2- The development of the halo was understood as changes in 
the occupation in the next shell or the sub-shell close to the continuum limit 



The central proton density is found to be decreasing near the neutron drip 



line, which is due to the proton-neutron interaction. However the diffuseness 

of the proton density does not change for the whole Na isotopes. 
PACS numbers : 21.10.Gv,21.60.-n,24.10.Cn,25.45.De 
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The recent developments in the accelerator technology and the detection techniques all 
around the world have changed the nuclear physics scenario. It is now possible to produce 
and study the nuclei far away from the stability line - so called "EXOTIC NUCLEI". 
Experiments of this kind have casted new light on nuclear structure and novel and entirely 
unexpected features appeared: e.g. the neutron halo in u Li [Q and neutron skin as the 
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rapid increase in the measured interaction cross- sect ions in the neutron-rich light nuclei. 
The extreme proton and neutron ratio of these nuclei and physics connected with these low 
density matter have attracted more and more attentions in nuclear physics as well as other 
fields such as astrophysics. 

With the exotic matter distribution near the drip line, a lot of questions are still open, 
e.g., the relation between the halo and the shell effect , the difference about the proton and 
neutron distribution on the stability line and away from the stability line. How is the halo 
formed ? Are there a rapid change from the normal nuclear density to the halo density or 
a gradual change in the particle number ? As the matter distribution is not measurable 
directly , series of experiment at different incident beam energy are necessary in order to 
determine the density distribution of both proton and neutron model-independently. Among 
all the experiments carried out so far, Na isotopes provide a good opportunity to study the 
density distributions over a wide range of neutron numbers 0. Although theoretically, 
lots of works have been reported with either the non-relativistic Hartree-Fock or relativistic 
mean field, but the pairing property is always neglected or simply treated by the BCS 
approximation, e.g., see [Q, ||. For the pairing correlation, the contribution from the 
continuum are essential for the description of the drip line nuclei. We report in this letter 
a systematic study of nuclear density distributions in Na isotopes within the Relativistic 
Mean Field Theory (RMF), with the pairing and the blocking effect for odd particle system 
properly described by Hartree-Bogoliubov theory in coordinate representation, and try to 
answer some of the general questions in these low-density nuclear matter. 

Definitely it is necessary to have a self-consistent theory to describe directly the cross 
section, so that the usual model dependent way for extracting the matter distribution could 
be avoided. As discussed in great detail in a recent review article [0 and in the references 
given therein, one has applied so far rather different techniques to describe halo phenom- 
ena in light nuclei, as for instance the exact solution of few-body equations treating inert 
sub-clusters as point particles, or the density dependent Hartree-Fock method in a localized 
mean field taking into account all the particles in a microscopic way, or a full shell model 
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diagonalizations based on oscillator spaces with two different oscillator parameters for the 
core- and the halo-particles. Recently, a fully self-consistent calculation within the Relativis- 
tic Hartree-Bogoliubov (RHB) theory in coordinate space for the description of the chain of 
Lithium isotopes ranging from 6 Li to n Li was reported [/]. It combines the advantages of a 
proper description of the spin-orbit term with those of full Hartree-Bogoliubov theory in the 
continuum, which allows in the canonical basis the scattering of Cooper pairs to low lying 
resonances in the continuum. Excellent agreement including binding energy, matter radius, 
and matter distribution with the experimental data was obtained without any modification 
of the neutron lpi/% level like the former works ( see |7j] and the references therein ). 

We first study the isospin dependence of the density distribution and the ground state 
properties of the Na nuclei within RHB. Then the cross sections based on Glauber model 
calculations with the density obtained from RHB are directly compared with the experi- 
mentally determined ones. As the theory used here is fully microscopic and parameter free, 
it gives consistent description of the proton and neutron distribution, and the development 
of proton and neutron halo or skin could be examined. 

The basic ansatz of the RMF theory is a Lagrangian density whereby nucleons are de- 
scribed as Dirac particles which interact via the exchange of various mesons and the photon. 
The mesons considered are the scalar sigma (a), vector omega (u) and iso- vector vector rho 
(p). The latter provides the necessary isospin asymmetry ||. The scalar sigma meson moves 
in a self-interacting field having cubic and quartic terms with strengths #2 and gs, respec- 
tively. The Lagrangian then consists of the free baryon and meson parts and the interaction 
part with minimal coupling, together with the nucleon mass M and m a (g a ), {gj) and 
m p (g p ) the masses (coupling constants) of the respective mesons: 

C = - M)V + §d M <70"<7 - U(a) - {tt^W 

+|m> M w" - \R^W V + \m 2 p p^ - \F pv F^ (1) 
-g^aip - g^ixpip - g p 4>pTip - eipfiip 

where a non-linear scalar self-interaction U(a) of the a meson has been taken into account 
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For the proper treatment of the pairing correlations and for correct description of the 
scattering of Cooper pairs into the continuum in a self-consistent way, one needs to extend 
the present relativistic mean- field theory to a continuum RHB theory ||10|| . Using Green's 
function techniques it has been shown in Ref. [O] how to derive the RHB equations from 
such a Lagrangian: 



(2) 



Ek are quasi-particle energies and the coefficients Uk(r) and Vk(r) are four-dimensional Dirac 
spinors. h is the usual Dirac Hamiltonian 
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containing the chemical potential A adjusted to the proper particle number and the meson 
fields a and to determined as usual in a self-consistent way from the Klein Gordon equations 
in no-sea-approximation. 

The pairing potential A in Eq. (||) is given by 

1 



A 



ab 



Vabcd K cd 



(4) 



cd 



It is obtained from the pairing tensor k = U*V T and the one- meson exchange interaction 
Vabcd m the pp-channel. As in Ref. M V^ cd in Eq. (|^) is the density dependent two-body 
force of zero range: 

V 



V( ri ,r 2 ) 



[l + P°)5{r x - r 2 ) (1 - p(r)/p ) 



(5) 



The RHB equations (0) for zero range pairing forces are a set of four coupled differential 
equations for the HB quasi-particle Dirac spinors U(r) and V(r). They are solved by the 
shooting method in a self-consistent way as [0. With a step size of 0.1 fm and using proper 
boundary conditions the above equations are solved in a spherical box of radius R = 25 fm. 



As shown in [12], for these relatively lighter nuclei R = 25 fm give quite accurate result. 



Since we use a pairing force of zero range @ we have to limit the number of continuum 
levels by a cut-off energy. For each spin-parity channel 20 radial wave functions are taken 
into account, which corresponds roughly to a cut-off energy of 120 MeV for R = 25 . After 
fixing the cut-off energy and the box radius R, the strength Vq of the pairing force (0) for 
both the neutrons and protons is fixed by a similar calculation of Gogny force D1S |]T3[ by 



reproducing the corresponding pairing energy — |TrA/ti as ref. 0. For p Q we use the nuclear 
matter density 0.152 fm -3 . The ground state > of the even particle system is defined 
as the vacuum with respect to the quasi-particle: (3 V \^ >= 0, \ty >= l\ u (3 U \— >, where 
| — > is the bare vacuum. For odd system, the ground state can be correspondingly written 
as: 1^ >n~ P^Hu^ii Pu\— >, where /i is the level which is blocked. The exchange of the 
quasiparticle creation operator M with the corresponding annihilation operator (3^ means 
the replacement of the column fi in the U and V matrices by the corresponding column in 



the matrices V*, U* fTS]] . 



A systematic set of calculations have been carried out for all the nuclei in Na isotopes 
with mass number A ranging from 17 to 45. We have employed in the calculations the 
non-linear Lagrangian parameter set NLSH which was widely used for the description of all 
the medium and heavy nuclei, particularly drip line nuclei ||TB| . 



The calculated binding energies Eb and the interaction cross sections with the Glauber 
Model are presented in Figure 1. The calculated binding energies Eb are in good agreement 



with the empirical values |17| for the radioactive isotopes, which are our current inter- 
ests. The resonance states of 17 Na and 18 Na ( with a positive Fermi energy ) are exactly 
reproduced. 19 Na is bound but unstable against the proton emission, reproducing the ex- 
perimental observation. The neutron drip-line nucleus has been predicted to be 45 Na in 
the present model. The difference between the calculations and the empirical values for the 
stable isotopes is from the deformation, which has been neglected here. 

To compare the cross section directly with experimental measured values, the densities 
Pn,p( r ) of the target 12 C and the Na isotopes obtained from RHB ( see Fig. 2 ) were used. 
The cross sections were calculated in Glauber model by using the free nucleon-nucleon cross 



section |14] for the proton and neutron respectively. The cross sections for reaction of Na 
isotopes at 950^4 MeV on 12 C have been compared with the experimental values |5) in the 
upper part of Fig.l. The agreement between the calculated results and measured ones are 
fine. The cross section below 22 Na changes only slightly with the neutron number, which 
means the proton density has played important role to remedy the contribution of less 
neutron. From 25 Na to the neutron drip line, a gradual increase of the cross section has 
been observed. After 32 Na, although no data exist yet, a relatively fast increase has been 
predicted. As we have seen here, the measured cross section shows similar behavior with 
that of RHB. 

As the density, which is obtained from a fully microscopic and parameter free model, 
is well supported by the experimental cross sections and binding energies, we proceed to 
examine the density distributions of the whole isotopes and study the relation between the 
development of halo and shell effect within the model. 

The density distributions for both the proton and the neutron are given in Fig. 2 . As 
seen in the upper part of Fig. 2 , the change of the neutron density is as follows: for the 
nucleus with less number of neutrons (N), the density at the center is low and it spreads 
only to some smaller distance. With the increase of N, the density near the center increases 
due to the occupation of the 2s\/2 level, so does the development of neutron radius. In the 
lower part of Fig. 2 , the proton density shows different behavior. The surface is more or 
less unchanged because of the Coulomb Barrier, but with the increase of N, the density of 
the center decreases due to the slight increasing at the tail. This is considered to be due 
to the attractive proton neutron interaction. But as the density must be multiplied by a 
factor Anr 2 before the integration in order to give the fixed Z, the big change in the center 
does not influence the outer part of the proton distribution very much. This result is fully 



consistent with the recent experiment on the charge- changing cross section JT8[ . 

The matter densities for the even N Na isotopes are given in Fig.3. the shortest tail 
in the total density occurs for 23_25 Na, the most stable ones. For either the proton or the 
neutron rich side, the tail density increases monotonically. The tail ( r > 10 fm ) of the 



proton rich nuclei is mainly due to the contribution of the proton and that of the neutron 
drip-line is mainly due to the contribution of the neutron. Compared with the neutron-rich 
isotopes, the proton distribution with less N has higher density at the center, lower density 
in the middle ( 2.5 < r < 4.5fm ), a larger tail in the outer ( r > 4.5 fm ) part. 

Next we will examine the density distribution for the neutron rich side. It is interesting 
to connect the matter distribution with the level distribution ( See Fig. 4 ): after 25 Na, as 
it is a sub-closure shell for the lcfs/2, then the neutrons are filled in the 2s\/2 and 1g?3/2- So 
the tail of the density for 27 Na is two order of magnitude larger than 25 Na at r = 10 fm, 
while the tail of the density from 27 Na to 31 Na is very close to each other. But as more 
neutrons are filled in, the added neutrons are filled in the next shell I/7/2, 2p 3 / 2 and 2pi/ 2 . 
So again two order of magnitude's increase has been seen from 31 Na to 33 Na, and then a 
gradual increase after 33 Na. So it becomes clear that the rapid increase in the cross section 
is connected with the filling of neutrons in the next shell or sub shell. In the inserts in 
Fig. 3 : the radius r at which p(r ) = 10~ 4 fm -3 is given as a function of the mass number 
to see the relation between the shell effect and matter distribution more clearly. It is very 
interesting to see a slight decrease of r from proton drip-line to 25 Na. The tail of 25 Na is 
the smallest. From 26 Na to 31 Na, one sees a almost constant r . After a jump from 31 Na to 
32 Na, a rapid increasing tendency appears again. 

In Figures 4 , the microscopic structure of the single particle energies in the canonical 



basis [[HJ is given. In the left panel of Fig. 4, the single particle levels in the canonical basis 
for the isotopes with an even neutron number are shown. Going from A = 19 to A = 45 
we observe a big gap above the N = 8, N = 20 major shell , and N = 14 sub-shell. The 
N = 28 shell for stable nuclei fails to appear, as the 2p 3 / 2 and 2pi/ 2 come so close to 1/7/2- 
When N > 20, the neutrons are filled to the levels in the continuum or weakly bound states 
in the order of I/7/2, 2p 3 / 2 , 2pi/ 2 , and 1/5/2- In the right part, the occupation probabilities 
in the canonical basis of all the neutron levels below E = 10 MeV have been given for 35 Na 
to show how the levels are filled in nuclei near the drip line. The importances of careful 
treatment of the pairing correlation, of treating properly the scattering of particle pairs to 



higher lying levels, are noted in the figure. 

Summarizing our investigations, the development of a proton skin as well as neutron 
skin has been systematically studied with a microscopic RHB model, where the pairing and 
blocking effect have been treated self-consistently. A systematic set of calculations for the 
ground state properties of nuclei in Na isotopes is presented using the RHB together with 
standard Glauber theory. The RHB equations are solved self-consistently in coordinate space 
so that the continuum and the pairing have been better treated. The calculated binding 
energies are in good agreement with the experimental values. A Glauber model calculation 
has been carried out with the density obtained from RHB. A good agreement has been 
obtained with the measured cross sections for 12 C as a target and a rapid increase of the 
cross sections has been predicted for neutron rich Na isotopes beyond 32 Na. The systematics 
of the proton and neutron distribution are presented. After systematic examination of the 
neutron, proton and matter distributions in the Na nuclei from the proton drip-line to the 
neutron drip-line, the connection between the tail part of the density and the shell structure 
has been found. The tail of the matter distribution is not so sensitive to how many particles 
are filled in a major shell. Instead it is very sensitive to whether this shell has occupation 
or not. The physics behind the skin and halo has been revealed as a spatial demonstration 
of shell effect: simply the extra neutrons are filled in the next shell and sub-shell. This is in 
agreement with the mechanism observed so far in the halo system but more general. As the 
I/7/2 is very close to the continuum, the N = 28 close shell for stable nuclei fails to appear 
due to the coming down of the 2p 3 / 2 and 2pi/ 2 levels. Another important conclusion here is 
that, contrary to the usual impression, the proton density distribution is less sensitive to the 
proton and neutron ratio. Instead it is almost unchanged from the proton drip-line to the 
neutron drip-line. Similar conclusion has been obtained recently by charge change reaction 
experimently JT3J. The influence of the deformation, which is neglected in the present 
investigation, is also interesting to us, more extensive study by extending the present study 
to deformed cases are in progress 
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Figure Captions 

Fig. 1 Upper part: The interaction cross sections 07 of A Na isotopes on a carbon target at 
950A MeV: the open circles are the result of RHB and the available experimental data 
( A = 20 — 23, 25 — 32 ) are given by solid circles with their error-bar. The dashed 
line is a simple extrapolation based on the RHB calculation for 28-31 Na. Lower part: 
Binding energies for Na isotopes, the convention is the same as the upper part, but 
the RHB result for particle unstable isotopes are indicated by triangle. 

Fig. 2 The neutron ( upper ) and proton ( lower ) density distributions in Na isotopes. 
The same figures but in logarithm scale are given as inserts to show the tail part of 
the density distribution. 

Fig. 3 The same as Fig. 2 but for matter density distribution. The upper part is given in 
logarithm scale ( the radius is labeled at the top of the figure ) and the radius ro at 
which p(r ) = 10~ 4 for different isotopes is given as inserts. 

Fig. 4 Left part: Single particle energies for neutrons in the canonical basis as a function 
of the mass number. The dashed line indicates the chemical potential. Right part: 
The occupation probabilities in the canonical basis for 35 Na. 
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